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We theoretically calculate the interaction-induced frictional Coulomb drag resistivity between two 
graphene monolayers as well as between two graphene bilayers, which are spatially separated by a 
distance "d". We show that the drag resistivity between graphene monolayers can be significantly 
affected by the intralayer momentum-relaxation mechanism. For energy independent intralayer 
scattering, the frictional drag induced by inter-layer electron-electron interaction goes asymptotically 
as pd ~ T^/n'^d'^ and po ~ /n^d^ in the high-density {kpd S> 1) and low-density [kpd <^ 1) 
limits, respectively. When long-range charge impurity scattering dominates within the layer, the 
monolayer drag resistivity behaves as po ~ /n^d"^ and hi{y/nd)/n for kpd ^ 1 and kpd <^ 1, 
respectively. The density dependence of the bilayer drag is calculated to be pn oc T"^ jrC" both in the 
large and small layer separation limit. In the large layer separation limit, the bilayer drag has a strong 
l/d* dependence on layer separation, whereas this goes to a weak logarithmic dependence in the 
strong inter-layer correlation limit of small layer separation. In addition to obtaining the asymptotic 
analytical formula for Coulomb drag in graphene, we provide numerical results for arbitrary values 
of density and layer separation interpolating smoothly between our asymptotic theoretical results. 
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I. INTRODUCTION 

Much attention has recently focused on multilayer 
systems in graphene, where carrier transport proper- 
ties may be strongly affected by interlayer interac- 
tion effectsi"— . In particular, temperature and den- 
sity dependent Coulomb drag properties have recently 
been studied in the spatially separated double layer 
graphene systems^. Frictional drag measurements of 
transresistivity in double layer systems have led to 
significant advances in our understanding of density 
and temperature dependence of electron-electron interac- 
tions in semiconductor-heterostructure-based parabolic 
2D systems2r— . While electron-electron interactions 
have indirect consequences (for example, through car- 
rier screening) for transport properties of a single isolated 
sheet of monolayer or bilayer graphene, the Coulomb drag 
effect provides an opportunity to directly measure the ef- 
fects of electron-electron interactions through a transport 
measurement, where momentum is transferred from one 
layer to the other layer due to inter-layer Coulomb scat- 
tering. In Coulomb drag measurements the role of elec- 
tron interaction effects can be controlled by varying, in 
a systematic manner, the electron density (n), the layer 
separation (d), and the temperature (T) since the in- 
terlayer electron-electron interaction obviously depends 
on n, d and T . (There is also a rather straightforward 
dependence of the drag on the background dielectric con- 
stant K, arising trivially through the interaction coupling 
constant or equivalently the graphene fine-structure con- 
stant, which we do not discuss explicitly.) 

In view of the considerable fundamental significance 
of the issues raised by the experimental observations^, 
we present in this paper a careful theoretical calcula- 
tion of frictional drag resistivity, p£){T), both between 
monolayer graphene (MLG) sheets and bilayer graphene 



(BLG) sheets, using Boltzmann transport equation. The 
current work is a generalization of the earlier theoretical 
work on graphene drag by Tse et ali^, where the drag is 
calculated within the canonical many-body Fermi liquid 
theory assuming an energy independent intralayer mo- 
mentum relaxation time. We would like to note that, al- 
though an energy independent relaxation time captures 
the effects of both short range and screened Coulomb 
impurities in 2D electron gas with parabolic dispersion, 
it does not correspond to any disorder model in MLG 
with its linear dispersion. In this paper, we generalize 
the formalism for calculating drag resistivity by includ- 
ing an arbitrary energy dependent intralayer scattering 
mechanism within the linear response Boltzmann equa- 
tion description. One new feature of our work is that 
we consider both MLG and BLG drag on equal theoreti- 
cal footing, providing detailed theoretical drag results for 
both systems. 

We find that, for MLG, the energy-dependent in- 
tralayer transport scattering time is the key to under- 
standing the density dependence of drag resistivity in- 
duced by inter-layer electron-electron interaction. In 
the presence of an energy independent scattering time, 
the low temperature drag goes asymptotically as pn ^ 
T'^/n'^d^ and po ~ T'^/n^d'^ in the high-density {kpd > 
1, where kp ^s. the Fermi wave vector) and low-density 
{kpd <^ 1) limits, respectively. These interlayer drag 
results for energy-independent intralayer momentum re- 
laxation, however, change qualitatively in the presence 
of an energy-dependent intralayer relaxation. For ex- 
ample, if the intralayer scattering is dominated by the 
charged impurities, which is believed to be the dominant 
intralayer momentum relaxation mechanism in most cur- 
rently available graphene samples on a substrate^^"^^ , the 
MLG drag resistivity behaves asymptotically as ^ 
T'^/n^d'^ and T'^\x\{^d)/n for kpd > 1 and kpd < 1, 
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respectively. In the low temperature regime {T /Tp ^ 1), 
the enhanced phase space for q = 2kp Coulomb backscat- 
tering leads to ~ lii(T) corrections to the usual de- 
pendence of the drag resistivity in ordinary 2D electron 
systems, i.e. po oc T^lnT. We find that, due to the 
chirality induced suppression of the q — 2kF backscat- 
tering in MLG, these ~ lii(T) corrections to the drag 
resistivity are absent in this system. 

We also investigate the Coulomb drag resistivity for 
two bilayer graphene sheets separated by a distance "d" . 
The low temperature behavior of BLG drag follows the 
usual dependence found in MLG systems. Both 
in the weakly correlated large layer separation limit, 
qTpd ^ 1, and the strongly correlated small separa- 
tion limit, qxpd 1 (qxp being the Thomas Fermi 
screening wavevector in BLG), the bilayer drag shows 
an inverse cubic dependence on the carrier density; i.e. 
Pd oc T^/n^. In the large layer separation limit, the BLG 
drag has a strong dependence on layer separation, 
whereas this goes to a weak logarithmic dependence in 
the strong inter-layer correlation (small layer separation) 
limit. The BLG drag, in contrast to MLG drag, does 
not manifest any qualitative dependence on the intralayer 
momentum relaxation mechanism with both short-range 
disorder scattering and long-range charged impurity scat- 
tering producing the same interlayer BLG drag. 

The paper is organized as follows. In Sec. II we present 
the general formula for drag resistivity in 2D materials 
from a Boltzmann transport approach, considering arbi- 
trary momentum dependent scattering times. The anal- 
ysis here is for arbitrary 2D chiral electron systems with 
two gapless bands, and hence is applicable to both MLG 
and BLG drag resistivity. In Sec. Ill, we study in de- 
tail, the temperature, density and layer separation depen- 
dence of the drag resistivity in MLG. We provide both 
analytic low temperature asymptotic forms and numeri- 
cal results for the MLG drag resistivity, using an energy 
independent as well as a linearly energy dependent scat- 
tering time. In Sec. IV, we study the drag resistivity of 
BLG, focusing both on analytic asymptotic forms and nu- 
merical results, within an energy independent scattering 
time approximation. Finally, we conclude our study in 
Sec. IV. with a summary of our work and a comparison 
of our results with other works in this field. 



II. DRAG RESISTIVITY IN 2D CHIRAL 
ELECTRON SYSTEMS 

In this section, we will derive the general formula for 
drag resistivity of chiral 2D electron-hole systems. We 
consider two layers (a and p) of the material, which are 
kept separated by an insulating barrier of thickness d, 
such that the carriers in the different layers are only cou- 
pled through the Coulomb interaction and there is no 



tunneling between the two layers. In drag experiments, 
an electric field is applied to layer a (the driven or ac- 
tive layer) and causes a current density ja, which induces 
the electric field Ep in layer p (the dragged or passive 
layer), where the current jpis set to zero. Then the drag 
resistivity is defined by = E^/j", where a is the di- 
rection along which the current ja flowsi^i^. 

We will consider each layer to be composed of a chiral 
two-band (electron and hole) material where the den- 
sity of the carriers can be changed independently. Let 
the dispersion of these bands be given by Sgk, and the 
corresponding two-component electron wavefunctions be 
given by ips'k, where s = ±1 denotes the band indices. The 
group velocity of the bands are given by v^k = Vk^sk- 

Let /^jj. denote the non-equilibrium distribution func- 
tion of the band s in the layer i, where {i = a,p). The 
Boltzmann equation for the distribution function is 

eE,.v,k^=i:k (1) 

where, I^k is the collision integral. Within a linearized 
relaxation-time approximation, 

BP (^ — P ) 
Jsk = /sk-eTkErv^k^— = Jsk + 'Psk tt; (^J 

CEsk 

where (/)*k = er^E^ • v^k, — tI the transport scat- 
tering time in the layer i, f^^^ = l/[e(^=''-A'i)/T _j_ the 
equilibrium Fermi distribution function in layer i, /ij the 
chemical potential in the layer i, T the temperature of 
the system. We have set fi. = 1 and the Boltzmann con- 
stant kg = \ throughout this paper. The current in layer 
i can then be written as 

= -ge 51 Vsk/^k = 5e X! '^kVskX^k (3) 

sk sk 

where 17 = 4 is the spin and valley degeneracy of the 
excitations. We will now focus on the form of the col- 
lision integral in the presence of the two layers. The 
collision integral in each layer has two contributions; (i) 
from impurity scattering in the same layer and (ii) from 
Coulomb scattering with electrons of the other layer; i.e. 
^Ik = ^t""^ and hence j. = + f\ The 

impurity scattering contribution to the current jlt"^^^ — 
cr'Ei, (T* being the usual conductivity of the layer i in 
absence of the second layer. From now on, we will focus 
on the Coulomb scattering contribution and drop the su- 
perscript (C) in the collision integral and the current. 

The electron-electron scattering between the layers 
is mediated through a dynamically screened inter-layer 
Coulomb interaction F(q, uj). We will discuss the precise 
form of the screened Coulomb interaction at the end of 
this section. The collision integral is given by 
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11^ = 2ngJ2J2l du;\V{q,u;)\'FC:;^Fi,r;_^ 

s'r,r' k'q 

X d{uj + Esk - es'-k+q)5{uj + Er'k'-q " £rk') (4) 

where r, r' = ±1, I is the layer other than i, and Fj^^^ = IV'I' k+qV''ikP is the wavefunction overlap between the bands. 

It is easy to veriiy that the collision integral vanishes if Using Eq. ([2]) , to linear order in the deviations from the 
we replace /' by the equilibrium approximation, = equilibrium distribution, the collision integral is given by 



(1 fsiJfs'k+qi^ /rk')/r'k'-q /sk(l /s'k+q)/rk' /r'k'-q) 



= ^ E E / du;\V{q,Uj)\'F^^^F^f_^f:^{l - /:,k+q)/'k'(l - /''k'-q)(Ck - '^^'k+q + 0'k' - 0''k'-q) 



s'r.r' k'q 

X 5{UJ + Esk - £s'k+q)<5(t^ + £r'k'-q ^ ^rk' ) (5) 

I 

Multiplying Eq. ([5]) by egr^v^k and summing over s and we get the current in layer z, 
k, the terms with 0* vanishes. Further, using the identity 

/sk(l - /s'k')'5('^ +£sk - £s'k') = ^^T^ — (6) 

1 — e T 

I 

J» = ^r;^ 2^ / du „ ^ 2^ ^k,qT-kVsk(/,k - /.'k')'5(w + e^k - e^'k+q) (7) 

^ J smh ss'k 

X E ^k^q^' ■ (^k'Vrk' - T-k'+qV,,k'+q)(./;k' " /r'k' + q)'5(^ + ^rk' - e^'k' + q) 
rr'k' 



Rearranging the terms with we finally obtain 



where 
Gd = 



AttT ^ J 



j? = <^D^?. (8) 



|^(q,c.)|2/m[x^(q,c^)]/m[x"(q,c<.)] 



sinh ^ 



(9) 



and the non-linear drag susceptibility is given by 

-Pkfq(/Jk - /l'k+q)('^kVsk - ■^k+qVs'k+q) 



t,(q,w) = eE 



ss'k 



+ esk - £s'k+q + ?'0+ 



Then the drag resistivity is given by 



PD 



(10) 



(11) 



We will now focus on the last piece of informa- 
tion required to calculate the drag resistivity of chi- 



ral electron- hole systems, the dynamically screened in- 
teraction, V{q,uj), as shown in Fig. [IJa). The intra- 
layer bare Coulomb interaction is given by 14a (q) = 
^ppil) = 27re^//cg, where k is the dielectric constant. 
The inter-layer bare Coulomb interaction is given by 
Vap{q) = Vpa{q) = Vaa{q)e~'^'^ , where d is the layer sepa- 
ration. Within random phase approximation (RPA), the 
dynamically screened inter-layer interaction is given by 
V{q,uj) = Vap{q)/e{q,LL!), with the dielectric function of 
coupled layer systems given byS. 

|e(q,^)| = [l-Vaa{q)Tla{q,Lo)][l-Vppiq)Ilp{q,Uj)] 

- Vap{q)Vpa{q)na{q,Uj)np{q,Uj), (12) 

where 11^ is the polarizability of the layer i. 



n,(q,w) = -ffE 



-^k,qUsk 



/s'k+q) 



iO + EsVl- £,s'k+c 



(13) 



Eq. ([5]) - Eq. thus completely defines the drag resis- 
tivity of chiral electron-hole systems in terms of the band 
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(a) 



(b) 




FIG. 1: (a) Screened interlayer Coulomb interaction in the 
RPA. The thin and thick lines are the bare and the screened 
interactions, respectively. The bare bubble represents the 
polarizability n(g,i:i;). (b) The leading order diagrams con- 
tributing to the drag conductivity. F indicates the non-linear 
susceptibility. 



dispersion and wavef unctions. In the next two sections, 
we will adapt these general equations to study the drag 
resistivity in monolayer and bilayer graphene. 

So far we have used the linearized Boltzmann equation 
to derive the drag conductivity which relates the induced 
electric field in one layer to the driving current in the 
other layer in a double layer system. Before moving on to 
the specific case of SLG and BLG, we would like to show 
that the linearized Boltzmann result captures the leading 
order result for drag conductivity within a diagrammatic 
expansion of the linear response Kubo formula. When 
an external field is applied to layer 1 and the induced 
current is measured in layer 2, the drag conductivity is 
given by the Kubo formulaiii 



1 

.,0 



4{t),J2{0) 



(14) 



where A is the area of the sample and Ji is the current 
operator in the z-th layer. 

The nonvanishing leading order diagrams correspond- 
ing to Eq. ([T4|) are given in Fig. [TJb). The two leading 
order diagrams can be written in a symmetric form 



q.iUJrr 



xV{q,U}jn + ^n)V{q,U}m), (15) 



where fin = l-KinT and ujm = 2TTimT are boson fre- 
quencies, V{q,uj) is the interlayer screened Coulomb in- 
teraction [Fig. [TJa)], and ri(g,w) is the three-point ver- 
tex diagrams (or the non-linear susceptibility) and given 
bj^-ia 

r(g, w) = T ^ Tr {G,„ G,„+<. J(g)a„+c.} 

+Tr {Ge„a„-<. J(g)a„-<.} , (16) 

where G^^ is the Green function, J(g) stands for the cur- 
rent and "Tr" the trace. To get the dc drag conductivity 



we need to perform an analytical continuation of exter- 
nal frequencies to a real value, iiln — >■ i^, and the limit 
51 — > should be taken. 

After summing over the boson frequencies, ujm, and 
performing an analytical continuation to a real value of 
n we have 

'^^ ri((j,w)r2(g,w) |V"(g,cj)|^ . 



— E 



IGnT'^J sinh^:^^ 

q 21 



(17) 

The nonlinear susceptibility with real frequencies is given 
by 



^^^'^^ " i y"t^e(tanh^-tanh^) 

X ^ Tr [{G; - Gt) G-+^J{p)Gt+^] 
p 

+ {{q,Lj)^{^q,~Uj)}, (18) 

where Gf = {e — H ±ij)^^ denotes the retarded (— ) and 
advanced (-f ) Green function for a given system with the 
Hamiltonian H, respectively. Here we use a damping 
constant 7 — l/2r to include the disorder scattering in 
the Green function. In the Boltzmann regime {lot 3> 1 
or kpl ^ 1, where kp is the Fermi wave vector and I — 
vpT Ya the mean free path), which corresponds to weak 
impurity scattering and the actual experimental regime 
of the high mobility graphene samples, we can treat the 
vertex correction in the current Jijq) within the impurity 
ladder approximationiSiiS. Then, the impurity-dressed 
current vertex becomes J = [rtr /T)dH / dk^ where Ttr is 
the transport time. By using the following equation 



G7+„(p)G++Jp) = 2rImG++Jp), 



(19) 



and expressing the matrix form of the Green function in 
the chiral basis, finally, we have the nonlinear suscepti- 
bility 

r(g,c^) = [J,{k)- Js'{k + q)]Fss'{k,k + q) 

ss' ,k 

xlm /'^fc ~ f^'k+q ^20) 

LU + €sk - es'k+q + «0+ 

Comparing the above equation with Eq. (|10|) we have 

T{q,Lo)^2lYax{q,u:). (21) 

Thus we recover the Boltzmann equation drag conductiv- 
ity result by employing the leading order diagrammatic 
expansion within the Kubo formalism. 



III. DRAG RESISTIVITY IN MONOLAYER 
GRAPHENE 

Monolayer graphene (MLG) is characterized by the 
presence of gapless linearly dispersing electron and hole 
bands with dispersion. 



Csk = vpk, 



(22) 
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where s = ±1 denotes the electron and hole bands and vp 
is the Fermi velocity of graphene. The linear dispersion 
of MLG implies that the group velocity of the bands are 
given by Wsk = "^Fk/lkj. Thus, contrary to usual semi- 
conductor systems, or the case of bilayer graphene to 
be treated later, the group velocity is constant in magni- 
tude and does not scale with the momentum of the band- 
states. As we will show, this leads to profound qualitative 
differences in the variation of the MLG drag resistivity 
with the carrier densities in the layers and the distance 
between the layers. The chirality of the electrons (holes) 
are encoded in the band wavefunctions 

^--7l( s ) (23) 

where 6'k is the azimuthal angle in the 2D k space. This 
leads to the wavefunction overlap factor F^"^ = (1/2) [1-1- 



ss' cos(6'k-Hq - fk)]- 



Before we go on to a detailed discussion of the non- 
linear drag susceptibility and the drag resistivity in MLG, 
let us first discuss the finite temperature polarizability in 
MLG which will control the finite temperature screening 
of the Coulomb potential. Although analytic expressions 
for graphene polarizability at T = has been worked out 
before^S, the finite temperature versions of the polariz- 
ability have not been calculated analytically. The full ex- 
pression of finite temperature polarizability is necessary 
to understand more precisely the temperature dependent 
drag including the plasmon enhancement effects^. Here 
we provide a finite temperature generalization of our ear- 
lier work on zero temperature graphene polarizability^^. 
The MLG polarizability of layer i is given by 



nj(g,w,; 



2-Bfj / TT 



■nv 



7 \8 - z2 

dx [/(x) + g{x)\ 



(24) 



{^^-^ + ZX + a;^^ sgn(a+) (^"l^ - zx + x'^^ sgn(a_) 

^ (z^-y^Y +(z2-y2)(^a: + x2) ^ (^'-fy +{z^-y2){-zx + X^) 



where x = k/kpi, y = q/kp, and z ^ u}/EF^, f{x) = [e-(^-^')/*' - gix) = [e'■'=+^'^)/^' + , a± = z^ - ± 2zx, 
with Epi and kpi being the Fermi energy and the Fermi wave- vector in layer i. Here, /i^ is the chemical potential in 
layer i in units of Epi (to be calculated self-consistcntly) an d tj =^ T/Epj. 



We now turn our attention to the non-linear drag sus- 
ceptibility and drag resistivity in MLG systems. The 
drag resistivity in MLG crucially depends on the vari- 
ation of the transport scattering time with energy (or 
equivalently momentum). In this paper, we will consider 
the drag resistivity of MLG using two different models of 
scattering time: (a) a momentum independent scatter- 



ing time Tk = r and (b) a scattering time scaling linearly 
with momentum (energy), Tk — To|k|, which results from 
the unusual screening of charge impurity potential in this 
linearly dispersive materiali*i^"— . 

We first consider the energy independent scattering 
time approximation. In this case, the non-linear drag 
susceptibility in MLG can be written as 



X,((?,w) = 4nvp ^ ( s 



ss'k 



k + q , k \ K^^qiflk+q - /I'k) 
|k + q\ |k| / a; - e^k+q + £s'k + «0+ 



(25) 



Within the energy independent scattering time approxi- 
mations, the intra-layer conductivities are given by cr^ = 
e^Ep,n/A. 

We focus on the analytic asymptotic behaviors of the 
drag resistivity in MLG at low temperatures, both in the 
large layer separation weak coupling limit (kpd^ 1) and 
in the small layer separation strong coupling \miit{kpd <C 
1). To calculate the asymptotic behavior of drag resis- 



tivity first we investigate the non-linear susceptibility for 
cj < vpq < Ep. Due to the phase-space restriction the 
most dominant contribution to the drag resistivity arises 
from Lu < Vpq at low temperatures. When we neglect the 
energy dependence in the transport times, i.e., Tk = r, 
then we obtain, for cj < vpq, 

X(g,a.)^^ — 7l-gV4fc|^- (26) 
ttEp Vpq V 
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With assumptions of a large inter-layer separation 
{kpd ^ 1, or qTpd ^ 1, with qte being the Thomas 
Fermi (TF) screening wave vector) and the random phase 
approximation (RPA) in which Ha is replaced by its value 
for the non-interacting electrons, we have the drag resis- 
tivity at high density and low temperature, 



In Fig. [2]^a) we show the calculated Coulomb drag as 
a function of temperature for an equal carrier density. 



Pd 



h 5!C(5) ikeTy 



e2 3-28 Ef^Ef, {qTFM<lTF,d){kF,dY{kFj? 

(27) 

where qxF = '^Tskp is the TF wave vector with the 
graphene fine structure constant = / kvf and C,{x) 
is the Riemann zeta function. This result shows that 
PD{n) cx n^^ and Pd{T) oc . For large layer separa- 
tion (i.e. kpd ^ 1) the back-scattering q « 2kF is sup- 
pressed due to the exponential dependence of the inter- 
layer Coulomb interaction Vi2{q) oc exTp{—qd)/q as well 
as the graphene chiral property. In this case the drag 
is dominated by small angle scattering and one expects 

For the strong interlayer correlation {kpd <^ 1) in the 
low-density or small-separation limit, the asymptotic be- 
havior behavior of drag resistivity becomes 



PD 



h 1 (fcBT)2 

6 EfiEf2 {kFid)(kF2d) 



(28) 



We have the same temperature dependence, Pd{T) ^ 
T^, but the density dependence becomes much weaker, 
Poin) ^ l/{nd)'^. At low densities (or strong interlayer 
correlation, kpd ^ 1) the exponent in the density de- 
pendent drag differs from -4. In an ordinary 2D systems, 
at low carrier densities and for closely spaced layers the 
backward scattering can be important since kpd ~ 1. 
In the low temperature range T/Tp <^ 1 the enhanced 
phase space for q = 2kF backward Coulomb scattering 
leads to ln(T) corrections to the usual dependence 
of the drag, i.e. pjj oc T^lnT. However, due to the 
suppression of the q = 2kp back-scattering due to the 
chirality of graphene there is no ln(r) correction in the 
drag resistivity of monolayer graphene. 

The drag resistivity within this approximation is plot- 
ted as a function of temperature and density in Fig. [2] 
for different layer separations. The parameters corre- 
sponding to the experimental setup of ref. |3 are used. 



ni 



"2 



10 cm~ , for two different layer separa- 
tions d = 50 A and d — 200 A. The overall temperature 
dependence of drag is close to the quadratic behavior, 
Pd oc T^. But we find a small corrections at low tem- 
peratures, especially at low values of kpd. In regular 2D 
systems there is a ln(T) corrections to the dependence 
of the drag. However, due to the suppression of the back- 
scattering in graphene such logarithmic correction does 
not show up in our numerical results except perhaps at 
extremely low temperatures. In Fig. [5] (b) the density 
dependent Coulomb drag is shown for different layer sep- 
arati ons. Our calculated Co ulom b drag resistivity follow 



a 10 




n (lO'^cm-: 



FIG. 2: The calculated drag resistivity by considering the en- 
ergy independent scattering approximation, (a) The temper- 
ature dependence of Coulomb drag for two different layer sep- 
arations d = 50A, d = lOOA and the equal electron densities, 
ni = n2 = lO^^cm"^. (b) The density dependent Coulomb 
drag for different layer separations d = 5, 10, 20 nm and at 
T = 200K. 



a n" dependence with a ~ — 2 at low carrier densities 
(or, kpd < 1), but as the density increases the exponent 
decrease to a ~ — 4. 

So far, we have considered the energy independent 
scattering time. However, it is known that the scatter- 
ing by the charged impurity disorder which inevitably 
exist in the graphene environment dominates and the 
scattering time due to the charged impurity is linearly 
proportional to the energy, n ^ e = rffc.— li^"— . In this 
approximation, the nonlinear drag susceptibility is given 

by 



X^ {q, w) = At^^vp ^ [s(k + q) - s'k] - 

'k 



k.qVJ sk+q 



esk^ 



(29) 



In this case, we will mainly focus on the low temperature the non- linear susceptibility for w < vpq is 
asymptotic for of the drag resistivity. For linearly energy 

dependent scattering time due to the charged impurities, ^ ^ 2r° hp u) 

TT vp Ep ^\ - q^/Akp 



Xdq,^) = -^ — — ^===^. (30) 
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FIG. 3: The calculated drag resistivity by considering the lin- 
early energy dependent scattering time, (a) The temperature 
dependence of Coulomb drag for two different layer separa- 
tions d = 50A, d = 200A and the equal electron densities, 
m = 712 = 10^^ cm~^ . (b) The density dependent Coulomb 
drag for different layer separations d = 5, 10, 20 nm and at 
T = 200if . 



Then the drag resistivity for a large inter-layer separation 
{kpd 1) is given by 



Pd 



h C(3) (kBT)^ 



23 EfiEf2 ikFtd){kF2d){qTFid){qTF2d)' 



(31) 

and for the strong interlayer correlation limit {kpd <C 1), 
we have 



PD = 



3 Ef, 



In 



2(9tFi + qTF2)d+l 



'^{iTFi + qTF2)d 



(32) 



Thus we see that when we consider the energy dependent 
intralayer scattering time we have very different asymp- 
totic behaviors compared to the results from energy in- 
dependent scattering time, i.e., we find for kpd 1, 
Pd ~ T'^/{n^d^) and for kpd <s: I, po ^ T^\n{^/^d). 
Thus in the presence of the strong charged impurity scat- 
tering the Coulomb drag resistivity follow a n" depen- 
dence with a < — 3 at high densities but as the den- 
sity decreases the exponent {a) increase. Based on our 
calculation we believe that the experimental departure 
from the n^^ behavior reported in Ref.— is essentially a 
manifestation of the fact that the asymptotic regime 
is hard to reach in low density electron systems where 
kpd ^ 1 hmit simply cannot be accessed. We predict a 
weak \n{n)/n density dependence in the low-density or 
small separation limit. 

In Fig. [Ha) we show the calculated Coulomb drag as 
a function of temperature for two different layer sepa- 
ration (a) d — 50 A (h) d = 100 A by considering the 
linearly energy dependent scattering time. The overall 
temperature dependence of drag increases quadratically 
and there is no logarithmic correction due to the sup- 
pression of the back-scattering. In Fig. Eljb) the density 
dependent Coulomb drag is shown for different layer sep- 
arations. The density dependent Coulomb drag follow 
a n" dependence with a ~ — 2. Based on our calcula- 
tion the consideration of the energy dependent scattering 
time is crucial to understand experiment measurements 
of Ref. IJ. 



IV. DRAG RESISTIVITY IN BILAYER 
GRAPHENE 

In this section we will study the drag resistance in a 
heterostructure made of two bilayer graphene (BLG) lay- 
ers separated by an insulating barrier and study its de- 
pendence on temperature, density and separation of the 
layers. Both BLG and MLG have chiral gapless electron 
and hole bands. However, BLG differs from MLG in two 
crucial aspects: (i) the bands have a parabolic disper- 
sion as opposed to linear dispersion in MLG and (ii) the 
chiral angle is double that in MLG leading to enhanced 
rather than suppressed backscattering between the quasi- 
particles. We will show in this section how these two 
differences lead to dramatic changes in the density and 
layer separation dependence of the BLG drag resistivity 
compared to MLG drag resistivity. 

BLG consists of an electron and a hole band with 
quadratic dispersion 



Esk — sk^ /2m 



(33) 



(s = ±1 for electron(hole) band) and the BLG mass is 
m — O.OSSrrie, where me is the electron mass. The group 
velocity in the two bands are v^k = sk/m, which scales 
linearly with momentum. The wavefunctions in these 
bands are given by 



V2\ ^ 



(34) 



which gives the overlap factor 



k,q 



(1/2)(1 



ss'q^ SVC? 0/|k -I- qp, where (f) is the angle between k and 

q- 

Without loss of generality, we will assume that the 
applied electric field, the induced electric field and the 
resultant currents are all in the x direction. We will also 
assume that the chemical potential in both layers is in the 
conduction (electron) band, i.e. they are electron doped. 
The chemical potential as a function of temperature is 
obtained by solving the integral equation 



g(feV2m-M.)/T ^ I 



+ 



1 



-(fcV2m+A'.)/T I 



1 



(35) 

where pi is the density of the layer i. It is an interesting 
feature of the BLG dispersion that in the non-interacting 
approximation, the chemical potential is independent of 
the temperature and is given by /i^ = Efi, where Efi is 
the Fermi energy of the electrons in layer i. 

We will first focus on the polarizability and hence the 
screened Coulomb interaction in BLG systems. The ana- 
lytic form for the BLG polarizability at zero temperature 
has been derived before^ ^, but we need to take into ac- 
count the temperature dependence of the screening to 
study the detailed behaviour of the drag resistivity with 
temperature. Here, we generalize our earlier results on 
BLG polarizability to finite temperature. The polariz- 
ability can be broken up into the intra-band [s = s' terms 
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in Eq. ([13])] and inter-band [s s' terms in Eq. (IT3l) ] lytically we obtain 
contributions. Working out the azimuthal integrals ana- 



gm 



nr™(q,^) = ^ dx 



if+x f-x) 



and 

nr*^'-(q,c.) 



Air Jq x(x^ + z) 

if+x ~ f-x) 



gm 



dx- 



47r Jq x{x'^ + z) 



\x^ -y"^]- {x^ + z) + sgniCi + 2xy) 



\x^ - + {x^ + z)- sgn{C2 + 2xy)^/C| - 4a;2y2 



+ iy,z ^ -y,-z) (36) 



+ {y,Z'^ -y,-z) (37) 



where y = q/kp,, x = k/kp-i., z = w/E'f,, Ci = ^-y^, C2 = z + ^x^+y^, and fsx^^l (e^''^" + 1) with = T/Ep, 

I 



We will now shift our attention to the non-linear drag 
susceptibility. For BLG systems, both charge impurity 
scattering and short range impurity scattering leads to a 
transport scattering time independent of momenta, and 



hence, we will only consider an energy independent scat- 
tering time, = T*, in this case. With this approxima- 
tion the nonlinear drag susceptibility is given by 



X. (g, c) = 1^ ^ [. (k + q) - .'k] ^^-.^ (38) 



ss'k 

where m is the BLG mass, m = 0.033me, where rrie is the free electron mass. 



The non-linear susceptibility Xi can be separated into tion. The azimuthal integrals can be done analytically 
an intra-band contribution and an inter-band contribu- and we get 



7- r x(intra) / \i 



T^kpiy cos ( 



dx — , ^ — ^ „— — (2 — > —z) 



l£_p x{x^ + z) ^4a;2y2 _ (^ - y2)2 



(39) 



and 



7- r x{inter) / \n 



T^kpi cose 
ny 



9(6) 



dx A^^,i^(2a;2 - z)V4x2y2_(22;2+2/2_;2)2 _ 



x(x2 — z) 



^ -z 



(40) 



where b = 2z — y^ and 0^ is the azimuthal angle related to the vector q. 



We note that the intra-band drag susceptibility is 
proportional to the intra-band polarizability of BLG 
(X™*™(q,a;) = (4r/m)qn™*™(q, cj)) only in the T = 
limit. This relation breaks down at finite temperatures. 
Finally, with a momentum independent scattering time, 
and a quadratic dispersion, the intra-layer conductivity 
can be written in the simple form cr* ~ riie^T^ /m. Note 
that the scattering times cancel in the expression for drag 
resistivity which is purely dominated by electron-electron 



interactions. In figure^] (a), we plot the temperature de- 
pendence of the BLG drag resistivity for two different 
layer separations with a common layer carrier density 
of ni = n2 — 10^2^771^2^ The drag shows a quadratic 
behaviour with logarithmic corrections at low tempera- 
tures. In figure Hfb), we plot the density dependence 
of the bilayer drag at T = IGOX for two different layer 
separations. 

We now focus on the low temperature asymptotic be- 



9 




FIG. 4: (a): The Coulomb drag resistivity in BLG as a function of temperature for two layer separations, d = 50A (thick black 
line) and d = 70 A (dashed red line). The density of the two layers ni = n2 = 10^^ cm~^. (b): The density dependence of 
bilayer Drag at a temperature T = lOOA". The thick black line is for d = 50A and the dashed red line is for d = 70A. We have 
used a dielectric constant n — 4 corresponding to boron nitride substrate 



haviour of the BLG drag resistivity. To obtain the 
leading temperature dependence oi po at low temper- 
atures, both X and 11 can be replaced by their intra- 
band contributions at T = in the limit of small q,uj, 
i.e. Im[x^{q,u!)] ~ {gm/2Tr)TiU}/kp, and the screened 
Coulomb interaction is replaced by its static value with 
Hi = 112 ~ 5"T-/(27r). The screening in BLG is con- 
trolled by the Thomas- Fermi wave- vector qtf = e^gm/n, 
where k is the background dielectric constant. For large 
layer separation, i.e. qTpd 3> 1, the screened interac- 
tion has the form V{q) ^ q/[DoqTF siiih{qd)] in which 
Dq = gm/2TT is the density of states of bilayer graphene 
at Fermi level, and the leading order drag-resistivity is 
given by 



PD 



1 7r2C(3) T2 



1 



e2 16 EfiEf2 {kFid){kF2d){qTFdf 



(41) 



We note that contrary to MLG, qxF in BLG is a con- 
stant and is independent of the density. Thus the large 
layer separation limit is not the same as the high den- 
sity limit (fcpd ^ 1). In the opposite limit of small 
layer separation {qTpd ^ 1) and strong interlayer cor- 
relations, the screened interaction takes the form V{q) ^ 
(qTF / Do)exp{—qd) / {q+2qTF) and the leading order drag 
resistivity is given by 



1 tt' 



T2 



PD 



Itf 



e2 24 EfiEf2 kFikF2 



[— In (7Tj;^d + 7 — 1 — ln4 + ...], 
(42) 

where 7 is the Euler constant. Thus the drag resistance 
~ both in the large and small layer separation limit. 
The density dependence of the coefficient of the term 
is l/(nin2)^^^ in both limits. In the large separation 
limit, the leading order term ~ whereas in the 

small separation limit, the drag resistivity shows a weak 
logarithmic dependence on the layer separation. We also 
note that, since backscattering is enhanced in BLG, there 
will be an additional In T correction to the formulae 
derived here. 



V. CONCLUSION 

In this paper, we have theoretically studied the fric- 
tional drag between two spatially separated MLG and 
BLG layers to lowest nonvanishing order in the screened 
interlayer electron-electron interaction using Boltzmann 
transport theory (which is equivalent to a leading-order 
diagrammatic perturbation theory). We find that the 
low temperature drag mostly shows a quadratic temper- 
ature dependence, both in MLG and BLG, regardless of 
the layer separation and density of carriers. However the 
density and layer separation dependence of the coefficient 
of the T2 term is very different for MLG and BLG. 

The density and layer separation dependence of low 
temperature MLG drag resistivity crucially depends on 
the variation of the intralayer momentum scattering time 
with energy. For energy independent intralayer scat- 
tering time (which does not correspond to any model 
of disorder in MLG, but captures the effects of both 
short range and screened Coulomb impurities in standard 
2DEG), the drag varies from po oc T^/{ndy for kpd ^ 1 
to pd oc T'^/(n'^dP) for kpd ^ 1. However, for energy 
dependent intralayer scattering times (corresponding to 
intralayer Coulomb scattering by random charged impu- 
rities in the environment) the power law of density depen- 
dence is significantly changed. Thus, an accurate mea- 
surement of interlayer MLG drag is in principle capable 
of distinguishing the main disorder scattering. In most 
currently available graphene samples the scattering due 
to charged impurity disorder dominates. In this case the 
intralayer scattering time depends linearly on the energy 
and the drag resistivity becomes pn oc hi{y/nd)/n for 
kpd <^ 1 and pn oc for kpd ^ 1. We note that 

the density dependence of drag resistivity is very sensi- 
tive to the experimental setup, so the density dependence 
does not have any universal power law behavior because 
one is never in any asymptotic regime with clear cut ana- 
lytical power law behavior. Experimental measurements^ 
therefore may not find any clear cut power law behavior 
in the density dependence of MLG drag since one is al- 
ways in the crossover regime. We also find that due to the 
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suppression of the q = 2kp back-scattering in graphene 
there are no ~ ^^(T) corrections to the low tempera- 
ture MLG drag resistivity. 

We would like to take this opportunity to compare 
this work, which is a more complete and updated ver- 
sion of ref. [22l . with recent works on low temperature 
MLG drag resistivit y -'^^i^^" —, which predict widely vary- 
ing density and layer separation dependence of the low 
temperature MLG drag resistivity. Within energy in- 
dependent scattering time approximation, Tse et ali^ 
studied the drag in MLG in the high density limit and 
obtained po ~ T'^/n^d'^. Katsnelson^- ^jgo obtained the 
same result in the high density limit {tpd 3> 1) and 
additionally considered the kpd <^ 1 limit of ref. |lj find- 
ing a pd ^ ln(nc?^)/n dependence in the low density 
limit. Narozhny^'^, on the other hand, found that the 
drag vanishes for linearly dispersive Dirac particles. In 
the current paper, within the energy independent scat- 
tering time approximation, we find that pu ^ T"^ /n^d^ in 
the high density limit and pn ^ T^/ri^d^ in the low den- 
sity limit. Narozhny2^ computes the drag to linear order 
in the inter-layer potential where it is known to vanish 
even for ordinary 2D electron systems. The drag resistiv- 
ity is at least quadratic in the inter-layer potential and 
its vanishing to linear order is a rather trivial result. The 
main difference between our present work and Ref. 13 and 
[25. arises because both these papers use a form of the non- 
linear drag susceptibility where the band group velocity 
scales linearly with the momentum k and results in a 
leading order q dependence of the vertex in the suscepti- 
bility. While this is true for quadratic band dispersions 
as in regular 2D and BLG systems, the group velocity in 
linearly dispersive MLG is constant in magnitude(wF), 
and combined with the overlap functions, gives rise to 
a leading order q^ momentum in the non-linear suscep- 
tibility in our case. This accounts for the discrepancy 
between our results and earlier results. 

Finally Peres et. aP* considers the drag resistivity of 
MLG within a linearly energy dependent scattering time 
approximation and obtains pu ~ T^/n'^d^ in the high 



density limit and pD ~ T^/n^d^ in the low density limit. 
These results match with our energy independent scatter- 
ing time approximation results, whereas within linearly 
energy dependent scattering time, we get pu ^ /n^d'^ 
in the high density limit and pD ^ ln(nd^)/rt in the 
low density limit. This discrepancy can also be under- 
stood in terms of the scaling of the current vertex factors. 
The vertex factor going into the drag susceptibility is 
fsk'Hc- Peres et. aU^ uses a drag susceptibility where the 
vertex factor in the susceptibility scales as q^ . Thus they 
recover our energy independent scattering time scaling. 
However, with a constant group velocity, as is the case for 
MLG, the energy dependent scattering time approxima- 
tion should lead to a leading order q dependence (coming 
from the momentum dependence of the scattering time) 
of the vertex and this accounts for the discrepancy be- 
tween our results and Peres et. al. The conceptual ele- 
ment of our work is the new role of group velocity in the 
nonlinear susceptibility defining the drag, which turns 
out to be crucial to the calculation of Coulomb drag in 
monolayer graphene (with its linear energy-momentum 
dispersion). We have also studied the drag resistivity in 
two spatially separated bilayer-graphene structures. We 
find that the drag resistivity shows a quadratic temper- 
ature dependence at low temperatures, as for monolayer 
graphene. The density dependence of the BLG drag is in- 
dependent of the layer separation with po T^/n^ both 
in the large and small layer separation limit. In the large 
layer separation limit (weak inter-layer correlation), the 
drag has a strong 1/d^ dependence on layer separation, 
whereas this goes to a weak logarithmic dependence in 
the strong inter-layer correlation limit. 
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